Magnetic order in the quasi-one-dimensional spin 1/2 chain, copper pyrazine dinitrate 
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We present the first evidence of magnetic order in the quasi-one-dimensional spin 1/2 molecu- 
lar chain compound, copper pyrazine dinitrate Cu(C4H4N2)(NC>3)2. Zero field muon-spin relax- 
ation measurements made at dilution refrigerator temperatures show oscillations in the measured 
asymmetry, characteristic of a quasistatic magnetic field at the muon sites. Our measurements 
provide convincing evidence for long range magnetic order below a temperature Tn = 107(1) mK. 
This leads to an estimate of the interchain coupling constant of \J'\/kB = 0.046 K and to a ratio 
\J'/J\ = 4.4 x 10~ 3 . 

PACS numbers: 75.10.Pq, 75.50.Ee, 76.75.+i 



Experimental studies of quantum many-body physics 
rely on the existence of simple and well characterised 
model systems. One such system is the linear chain 
Heisenberg antiferromagnet (LCHAFM) which is de- 
scribed by the Hamiltonian 



H = ^ [- JSiSi+i + ,9^bB • Si], 



(1) 



where J is the nearest neighbour exchange constant along 
the chain and B is an applied magnetic field. For B = 
0, the system will not show long-range magnetic order 
(LRO) for T > 0, but will display excitations in the form 
of S = 1/2 spinons which form a gapless, two-particle 
continuum 0, The magnetic field B can tune the 
system through a quantum critical point (QCP) Q, (i.e. 
a zero temperature phase transition) at |B| = \2J/gfiB\ 
where the spin state becomes fully polarized 0,0. 

Experimental realisations of the LCHAFM have 
been made in several copper chain compounds, the 
most successful of which is copper pyrazine dinitrate 
Cu(C 4 H 4 N2)(N03)2 or CuPzN. This material has a rela- 
tively low value of J(= 10.3(1) K, determined from high 
field magnetization and specific heat measurements |6|]), 
ensuring the QCP is at an experimentally accessible mag- 
netic field, but also a small ratio of intrachain to inter- 
chain exchange constants (estimated as \J'/J\ < 1CP 4 
0), showing that the chains are well isolated. This has 
allowed the magnetic field dependence of the excitations 
in the two-spinon continuum to be measured and led to 
the identification of an extended critical state [§|, con- 
firming long-stan ding pre dictions of theoretical and nu- 
merical studies [UEiE]]- Three-dimensional (3D) LRO 
has not been detected in CuPzN down to 70 mK on the 
basis of specific heat measurements .12]- 

However, as we 

shall demonstrate, implanted positive muons are a much 
more sensitive probe of 3D LRO than specific heat for the 
case of very anisotropic magnetic chains. In this letter we 



show for the first time that CuPzN does indeed undergo 
a transition into a state of 3D LRO and provide an esti- 
mate of \ J'/J\ which is larger than previously thought. 

Zero-field muon-spin relaxation (ZF /i + SR) measure- 
ments were made on a powder sample of CuPzN using the 
LTF instrument at the Swiss Muon Source (S/xS), Paul 
Scherrer Institut, Villigen, Switzerland and the MuSR in- 
strument at the ISIS facility, Rutherford Appleton Lab- 
oratory, UK. In a /i + SR experiment [f3T | spin polarised 
muons are implanted into the sample. The quantity of 
interest is the asymmetry A(t), which is proportional to 
the spin polarization of the muon ensemble. 

Example ZF /i + SR spectra measured using a dilution 
refrigerator at S/xS are shown in FigQ] In spectra mea- 
sured below a temperature T = 110 mK oscillations in 
the asymmetry are observed (Fig^a)). These oscilla- 
tions are characteristic of a quasi-static local magnetic 
field at the muon stopping site, which causes a coherent 
precession of the spins of those muons with a component 
of their spin polarization perpendicular to this local field 
(expected to be 2/3 of the total polarization). The fre- 
quency of the oscillations is given by vi = j^Bi/2n, where 
7 M is the muon gyromagnetic ratio (= 2ir x 135.5 MHz 
T _1 ) and Bi is the average magnitude of the local mag- 
netic field at the iih muon site. Any fluctuation in mag- 
nitude of these local fields will result in a relaxation of the 
oscillating signal, described by relaxation rates A^. The 
presence of oscillations at low temperatures in CuPzN 
provides convincing evidence that this material is mag- 
netically ordered below 110 mK. 

Two separate frequencies were identified in the low 
temperature spectra, corresponding to two magnetically 
inequivalent muon stopping sites in the material. We 
note here that two muon sites have been observed in 
other copper pyrazine-based compounds 0|. The pre- 
cession frequencies, which are proportional to the inter- 
nal magnetic field as experienced by the muon, may be 
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FIG. 1: ZF /i + SR spectra for CuPzN at three temperatures, (a) T = 36 mK: Oscillations are observed in the asymmetry, 
characteristic of long-range magnetic order. A fit is shown to Eq.(|2J. (b) T = 120 mK: the spectra are described by two 
components. A fast relaxing signal, probably due to the existence of a radical state, and a larger component, seen to be 
exponential in form. A fit is shown to Eq.lPJ. (c) T = 405 mK: The fast relaxing component is still evident, although the 
larger component now shows a Gaussian form. A fit is shown to Eq.Q. 



viewed as an effective order parameter for these systems. 
In order to extract the temperature dependence of the 
frequencies, the low temperature data were fitted to the 
functional form 

A(t) = Ax exp(-Ai^) cos(27rz/i£ + 4>{) 
+A 2 cxp(-A2^) cos(27w 2 i + 02) 
+A 3 exp(-At) + A hg exp(-Abg^), (2) 

where A3 exp(— At) accounts for the contribution from 
those muons with a spin component parallel to the lo- 
cal magnetic field. The term A\> g exp(— Ab g i) reflects the 
small signal from those muons which stop in the silver 
sample holder or cryostat tail (with Ab g <C A). 

Across the measured temperature range, the two fre- 
quencies were found to be in the proportions v\ : V2 = 
1 : 0.65 while the relaxation rates Xi were found to be in 
the ratio Ai : A2 = 1 : 0.75. These quantities were fixed 
in these ratios during the fitting procedure. The ampli- 
tudes and phases best describing the data are given in 
[l5j |. We note that the ratio of longitudinal to trans- 
verse components is in excess of the expected value of 
^3/(^1 + A 2 ) = 1/2 (see below). Nonzero phases were 
also required to fit the data, as observed in previous stud- 
ies of Cu-based chain compounds of this sort [14), The 
magnitudes of the frequencies were fitted as a function of 
temperature with other parameters in Eq.J2J fixed 
The resulting temperature evolution of the precession fre- 
quencies is shown in Figl^a). 

From fits of v { to the form Vi(T) = z^(0)(l - T/T^f 
close to Tn, we estimate Tn = 107(1) mK, leading to a 
value of the ratio \k-QT-^/ J\ = 0.0103(1). The parameter 
/3 is highly sensitive to the chosen value of Tn due to 
the difficulty of obtaining data sufficiently close to the 
transition. Using Tn = 107(1) mK, we estimate f3 = 



0.18(5), which is consistent with low dimensional (i.e. 2D 
Ising or 2D XY) behaviour. Fits to v t {T) = ^(0)(1 - 
(T/Tn)")' 3 with a = 3 yield 2/1 (0) = 1.922(4) MHz and 
^2(0) = 1.257(3) MHz corresponding to local magnetic 
fields at the two muon sites of B\ — 141.8(3) G and 
B 2 = 92.7(2) G. 

Further evidence for a magnetic transition is seen in 
the transverse relaxation rates Ai, which increase as the 
magnetic transition is approached from below (Fig|^b)), 
as would be expected for critical slowing of fluctuations 
near a magnetic transition. The longitudinal relaxation 
rate A decreases as the transition is approached from be- 
low (Figl^c)), suggesting that this quantity is dominated 
by the magnitude of the local magnetic field at the muon 
site. 

At temperatures just above Tn the measured spectra 
are best described by two components (Fig^b)). One is 
a fast relaxing contribution, with an amplitude ~ 5 %, 
likely to arise due to a fraction of the muons which ini- 
tially form radical states. This component will also exist 
in the magnetically ordered regime, but is not resolvable 
from the longitudinal relaxation and partly explains why 
the value of A3 is greater than that expected. The sec- 
ond component observed in this regime is an exponential 
functional form, characteristic of fast dynamics in the 
local magnetic field at the muon sites At higher 

temperatures still (T > 400 mK), this second compo- 
nent in the spectra appears Gaussian in form (Fig^c)), 
characteristic of the early time behaviour of a Kubo Toy- 
abe function due to the magnetic field distribution 
of a static, disordered system, with slow dynamics pre- 
venting the characteristic recovery of the polarization at 
late times. It is likely that at these elevated temper- 
atures, the electronic moments are fluctuating too fast 
compared to the muon time scale, and are therefore mo- 



3 




20 40 60 80 100 120 
T (mK) 



FIG. 2: Temperature evolution of the variable parameters in 
Eq.@ for T < 140 mK. (a) Muon precession frequencies Vi 
with fits shown to the expression m{T) = ^(0)(1- (T/Tn)")' 3 
(see main text), (b) Transverse relaxation rates \i are seen to 
increase as the magnetic transition is approached from below, 
(c) The longitudinal relaxation rate A, which decreases as the 
transition is approached from below. 
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FIG. 3: (a) Results of fitting data for T > 110 mK to Eq.©. 
(Dashed lines show T = Tn.) Inset: The parameter 8 chang- 
ing from 8 — 1 (exponential relaxation) to 8 — 2 (Gaussian 
relaxation) as the dynamics move out of the muon time win- 
dow with increasing temperature. Main panel: The relaxation 
rate As increases as the transition is approached from above, 
(b) LF /i + SR spectra made at T — 2 K in several applied 
magnetic fields Bl. Fits are shown to the LFKT function 
with A fixed at 0.271 MHz for all fields. 



tionally narrowed from the spectra. Instead the muon 
ensemble is depolarized by the magnetic field distribu- 
tion at the muon site due to the quasi-static nuclear 
moments of the surrounding nuclei. This interpretation 
is supported by the observation that applied magnetic 
fields of magnitude > 100 G repolarize the relaxation, as 
expected for a static distribution (repolarization implies 
Hmt_, 00 -4(i) = A(0)). The measured data for T > T N 
were therefore fitted to the function 

A(t) = A 4 exp(-A 4 i) + A 5 exp(-X 5 t) 5 + A hg (3) 

where the stretching parameter S describes the passage 
from exponential to Gaussian behaviour. The results of 
this analysis is shown in FigEJJa). The fact that the re- 
laxation rate A5 increases sharply as the proposed transi- 
tion temperature is approached from above lends further 
weight to the case for the existence of a magnetic transi- 
tion. 

In order to further investigate the high temperature 
behaviour of the material measurements were made in 
an applied longitudinal magnetic field (LF) Bl using a 
4 He cryostat at the ISIS facility. The pulsed structure of 
the muon source at ISIS has the advantage of allowing 



the observation of the muon polarization at longer times 
than at a continuous source such as SfiS. This enables 
us to better characterise the data with the Kubo Toy- 
abe function generalised for applied LFs . The loss of 
time resolution at a pulsed source causes the fast relaxing 
component described above to be missing from the sig- 
nal entirely. The results of these measurements, made at 
T = 2 K, are shown in Fig[3Jb). The spectra are found to 
be approximately temperature independent in the range 
2 < T < 50 K. The application of a field repolarizes the 
slowly relaxing component of asymmetry for B < 100 G 
consistent with the interpretation that the relaxation at 
these temperatures is due to the disordered, quasistatic 
nuclear magnetic moments. The measured ZF and LF 
spectra for both materials at 2 K are well described by 
the longitudinal field Kubo Toyabe (LFKT) function de- 
scribing the relaxation due to an ensemble of disordered 
magnetic moments, although the LFKT function must 
be multiplied by an exponential function with constant 
relaxation rate to account for the weak relaxation due 
to dynamic fluctuations [l7|. This treatment has only 
one free parameter A = (Bf oc ), where (Bf oc ) is the 
second moment of the average local magnetic field dis- 
tribution at the muon sites in ZF. This analysis yields a 
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TABLE I: Parameters for copper-chain c omp ounds Sr 2 Cu03 
[UEI, KCuF 3 |p, CPC(=Cu(py) 2 Cl2) Ea and CuPzN (0 
and this work) with values of \J' / J\ estimated from Eq.JlJ. 





\J\/ka (K) 


Tn (K) 


|fc B T N /J| 


\J'/J\ 




Sr 2 Cu0 3 


2200 


5.4 


2.5 x 10~ a 


9.3 x 10" 


-4 


KCuF 3 


406 


39 


9.6 x 10" 2 


5.2 x 10" 


-2 


CPC 


26 


1.1 


4.2 x 10~ 2 


2.1 x 10" 


-2 


CuPzN 


10.3 


0.107 


1.0 x 10" 2 


4.4 x 10" 


-3 



value of A = 0.271 MHz for CuPzN. 

An efficacious method for estimating the interchain 
coupling J' in LCHAFM's has recently been developed, 
based on a modified random phase approximation, mod- 
elled with classical and quantum Monte Carlo simula- 
tions |is||. This approach leads to an empirical formula 
for the interchain coupling, given by 



\J'\/k B = 



In 



(4) 



with A = 2.6 and c = 0.233. Substituting our value of 
T N = 107(1) mK and using J/k B = -10.3(1) K, yields 
\J'\/k B ~ 0.046 K and a ratio \J'/J\ = 4.4 x 10~ 3 . 
The parameters describing a number of copper-chain 
compounds are presented in table [I] Comparing val- 
ues of the ratio \J'/J\ shows that CuPzN is approxi- 
mately five times better isolated than copper bispyridinc 
dichloride (Cu(py) 2 Cl 2 or CPC) USUI, while I J '/ J \ for 
CuPzN is around one quarter of that for the best isolated 
LCHAFM, Sr 2 Cu0 3 (although this material has a QCP 
at a field around 200 times larger than that for CuPzN) 

num. 

The failure of specific heat studies |l2j to detect an 
anomaly due to the onset of 3D LRO in CuPzN down to 
70 mK, may be attributed to the value of the ratio \ J' / J\. 
Investigations using a stochastic series quantum Monte 
Carlo method J2f| show that the height of the expected 
peak decreases with decreasing \J'\, becoming nearly lin- 
ear below \J'/J\ = 0.2. For smaller value of the ratio 
\J'/J\, the entropy change at the critical temperature is 
so small that it makes an undetectable contribution to 
the measured specific heat. 

In conclusion we demonstrated the existence of long 
range magnetic order in the LCHAFM material CuPzN 
with Tn = 107(1) mK. This leads to an estimate of the in- 
terchain coupling of | J'l/fcs — 0.046 K. For many ID sys- 
tems, small values of \ J'/J\ prevent the detection of 3D 
LRO with specific heat, and reduced ordered moments 
make the use of other techniques problematical (coupled 



chain mean-field theory |2J] then predicts an ordered mo- 
ment size of ~ 0.05/xb for CuPzN). Thus we expect that 
muons may provide a particularly valuable tool for future 
magnetic measurements of the key parameters of highly 
anisotropic low-dimensional systems. 
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